Ultra-long range correlations of the dynamics of jammed soft matter 
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(Dated: June 30, 2010) 

We use Photon Correlation Imaging, a recently introduced space-resolved dynamic light scattering 
method, to investigate the spatial correlation of the dynamics of a variety of jammed and glassy 
soft materials. Strikingly, we find that in deeply jammed soft materials spatial correlations of 
the dynamics are quite generally ultra-long ranged, extending up to the system size, orders of 
magnitude larger than any relevant structural length scale, such as the particle size, or the mesh 
size for colloidal gel systems. This has to be contrasted with the case of molecular, colloidal and 
granular "supercooled" fluids, where spatial correlations of the dynamics extend over a few particles 
at most. Our findings suggest that ultra long range spatial correlations in the dynamics of a system 
are directly related to the origin of elasticity. While solid-like systems with entropic elasticity exhibit 
very moderate correlations, systems with enthalpic elasticity exhibit ultra-long range correlations 
due to the effective transmission of strains throughout the contact network. 

PACS numbers: 



I. INTRODUCTION 

Soft materials such as colloidal suspensions, emulsions 
and surfactant phases typically exhibit increasingly slow 
relaxation dynamics as a result of particle crowding or 
because of strong interactions, either attractive or re- 
pulsive jlj. Work in the past years has shown that a 
feature shared by most of these systems is the heteroge- 
nous character of their slow dynamics, which result from 
rearrangements that are localized in space and intermit- 
tent in time [l|-[lo| , in analogy with molecular glass form- 
ers [ill, [l3 driven athermal grains and foams 

For colloidal hard spheres, probably the most stud- 
ied soft matter model system exhibiting a glass tran- 
sition [Tsl, the size ^ of dynamical clusters undergo- 
ing correlated rearrangements has been shown to grow 
with volume fraction on approaching the glass transi- 
tion [1, 0, [1^ [lO]. For supercooled samples whose dy- 
namics are stationary, the growth is however modest, the 
largest reported values of ^ being of the order of a few 
particle sizes d, H, [l^ HO] • Similar results have been re- 
ported for weakly attractive systems 0, [1| ■ This agrees 
with numerical and experimental findings for molecular 
glass formers [HI [H, HI [13 and grains [ll4il|. Although 
the debate is still very active on whether or not the glass 
and the jamming transitions coincide in thermal hard 
spheres [23l - [28| , and thus on how far the regime where ^ 
grows may extend, it is unlikely that significantly larger 



'Electronic address: ' lucacip® Icvn. univ-montp2 . fr | 



correlation lengths may be measured in the supercooled 
regime. 

Recent experiments using time-resolved light scatter- 
ing methods suggest that this scenario could be very 
different for a wide range of soft materials quenched 
in a nearly-arrested, out-of-equilibrium state, which we 
shall refer to (somehow loosely) as jammed soft sys- 
tems. In these experiments, very large temporal fluc- 
tuations of the intensity autocorrelation function were 
observed 0, [13, [1^ , suggesting that the size of regions 
undergoing correlated rearrangements may be a sizeable 
fraction of the scattering volume, i.e. that they may ex- 
tend over macroscopic distances. Indeed, direct measure- 
ments of ^ in a strongly attractive colloidal gel ^3Qi] and in 
concentrated soft particles [3l| have shown that in these 
systems the correlation length of the dynamics is limited 
essentially only by the system size. 

In this paper, we present data on a variety of soft 
jammed systems, showing that extremely long ranged 
correlations of the dynamics in jammed systems are the 
rule rather than the exception. The systems investigated 
include hard and soft spheres, colloidal gels made of at- 
tractive particles, biomimetic protein films, a concen- 
trated surfactant solution ("onion" phase) and Laponite 
suspensions: with the exception of the supercooled hard 
spheres, ^ always exceeds 1 millimeter, much larger than 
any structural length scale. We discuss the role of both 
the strength and the microscopic origin of the elastic- 
ity in shaping spatial correlations of the dynamics and 
compare our results to the behavior of jammed materials 
under shear l32H35i. 
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II. MATERIALS AND METHODS 



A. Time resolved correlation and photon 
correlation imaging 

Although dynamic light scattering (DLS) is now a pop- 
ular and well established technique to probe the dynam- 
ics of soft materials , its usefulness to measure dynam- 
ical heterogeneity has been limited until recently by the 
spatial and temporal averages usually involved in DLS 
measurements. In this section, we briefly recall the main 
features of a series of DLS time- and space-resolved meth- 
ods that we have introduced in the past years. More 
details can be found in Refs. [4, 29, 30, 37, 38]. 

Time resolved correlation (TRC). Valuable informa- 
tion on spatial correlations of the dynamics can be ob- 
tained in time-resolved measurements, even if they lack 
spatial resolution. Intuitively, this is due to the fact that 
temporal fluctuations of the dynamics are enhanced if 
the probed sample volume contains a limited number 
of statistically independent regions, as is the case when 
the dynamics are correlated over large distances [T6| . 
More formally, a dynamical susceptibility can be intro- 
duced [3£i] , which quantifies temporal fluctuations of the 
dynamics, and which is proportional to the volume inte- 
gral of the spatial correlation of the dynamics, G4 j39l . l40| . 
In our DLS experiments, we achieve temporal resolution 
by using a CCD camera as a detector and by averag- 
ing the intensity correlation function over pixels rather 
than over time. The CCD is placed in the far field, so 
that each pixel is illuminated by light issued from the 
whole scattering volume at a well defined scattering an- 
gle 9. The dynamics are probed on a length scale £ ^ 1/q, 
where q = 47rnA~^ sin(6'/2) is the scattering vector and 
n and A are the solvent refractive index and the in- vacuo 
wave length of the laser that illuminates the sample, re- 
spectively. We measure a time-resolved intensity correla- 
tion function (proportional to the square of the dynamic 
structure factor), defined as 



sample 



Cl{t,T) 



{Ip{t)Ip{t + T)),^ 



(1) 



which we refer to as the degree of correlation. In Eq. ([TJ , 
Ip indicates the scattered intensity measured by the p- 
th pixel and the averages are taken over the whole CCD 
detector. The average dynamics and dynamical hetero- 
geneity are quantified by 

g2{T)-l=<cj{t,T) >t (2) 
X(t) = var[ci{t, r)] =< c/(i, r)^ >t - < ci{t, r) >2 , (3) 

where averages are taken over time and where g2 — 1 
is the usual intensity correlation function measured in 
DLS, while x is the equivalent for DLS of the dynamical 
susceptibility Xi used in numerical works. Note however 
that, contrary to X4, x is not normalized with respect 
to the number of particles in the scattering volume, a 
quantity not always easily accessible experimentally. 
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FIG. 1: Scheme of a PCI setup in two different configurations: 
a) 90° scattering angle; b) backscattering. From Ref. [s^ with 
permission. 



Photon correlation imaging (PCI). Spatial resolution 
may be achieved by modifying the collection optics so 
as to form an image of the scattering volume onto the 
CCD detector. This method, which we have termed pho- 
ton correlation imaging (PCI), has been described for a 
low-q setup in Ref. 30] . In Fig. [T] we show its imple- 
mentation for 6 — 90 (a) and 180 (b) degrees, as in the 
experiments reported here. In both cases, an image of the 
illuminated sample is formed onto the CCD with magni- 
fication M ~ 1, using light scattered within a small solid 
angle centered around a well defined scattering angle. 
The resulting image has a speckled appearance similar 
to that of conventional far-field scattering experiments. 
In practice, one adjusts the diaphragm aperture so as to 
roughly match the speckle size to the pixel size, which 
is typically of the order of 10 //m. Under these condi- 
tions, the speckle size is larger than the typical size of 
the scatterers (< 1 /im), which thus can not be resolved 
individually. Information on the local dynamics, how- 
ever, can still be obtained, because the relative motion 
of the scatterers results in fluctuations of the intensity 
of each speckle, as in regular DLS. In contrast to con- 
ventional DLS, here each speckle is illuminated by light 
issued from a well defined, small region of the sample, so 
that spatially-resolved dynamics can be measured. 

In practice, we divide the CCD images in regions of 
interest (ROIs), each corresponding to a small volume 
in the sample, and apply the TRC method separately to 
each ROI. We compute a space- and time-resolved degree 
of correlation, defined as 



c/(i,r,r) 



{Ip{t)Ipit + T))^ 



(4) 



where the average is now over all pixels within a ROI 
corresponding to a small volume centered around r. The 
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ROIs must contain at least about 100 speckles for the sta- 
tistical noise on the local degree of correlation to be ac- 
ceptable: this limits the resolution of the method, which 
should be regarded as a coarse-grained technique. The 
main advantage of PCI is to decouple the size of the 
field of view (dictated by the magnification M) from the 
length scale £ over which the dynamics are probed (dic- 
tated by q and thus the scattering angle). Thus, very 
restrained motion can be measured with a large field of 
view. This is not the case for conventional imaging meth- 
ods, where ^ is a fixed fraction of the field of view. 

The spatial correlation of the dynamics may be quan- 
tified by comparing the temporal evolution of the local 
degree of correlation calculated, for a given lag, for dif- 
ferent locations. Traces corresponding to pairs of re- 
gions with correlated dynamics will exhibit similar fiuc- 
tuations, while uncorrelated ROIs yield independent fiuc- 
tuations of c/ . An example that will be discussed towards 
the end of this paper is given in Fig. [HI We define 



hard spheres, an "onion" gel, and a Laponite suspen- 
sion. For the sake of comparison, we recall also results 
obtained from concentrated soft spheres and a colloidal 



G4(r, Ar) 



{6ci{t,T, r)Sci(t,T,r + Ar))f 
^ var [d (t, T, r)] var [a (t , t, r + Ar)] 



(5) 

where Sci = cj— < ci >t are the temporal fiuctuations 
of the local dynamics. This is the analogous, albeit at 
a coarse grained level, of the spatial correlation of the 
dynamics calculated in numerical and experimental work 
where particle trajectories are accessible [13, 19, 32, 40]. 
In most cases, the dynamics are isotropic and we average 
G4 over all orientations of Ar. 

An important point concerns the normalization of 
Eq. ([5]): we recall that c/ contains a noise contribution 
due to the statistical noise associated with the finite num- 
ber of pixels processed for each ROI [H]. In PCI exper- 
iments, this contribution may be quite large, since one 
typically tries to reduce the size of the ROIs as much as 
possible in order to achieve a better spatial resolution. 
Distinct ROIs have uncorrelated noise; as a result, the 
numerator of the r.h.s. of Eq. ([S]) is noise-free (except for 
Ar — 0), while the noise contributes to the denomina- 
tor by increasing its value. Thus, at spatial lags > G4 
is depressed because of the noise contribution, an effect 
that depends on the size of the ROIs used for the anal- 
ysis. For systems that exhibit a finite spatial correlation 
of the dynamics at Ar > 0, we introduce a normalized 
spatial correlation of the dynamics defined as 



G4(T,Ar) =6(T)G4(T,Ar), 



(6) 



where 6 is a time-lag dependent coefficient chosen so that 
G4 ^ 1 for Ar ^ 0. 



B. Experimental systems 

We report below PCI measurements of heterogeneous 
slow dynamics on a variety of systems: "artificial skin" 
(Vitro-corneum® by IMS Inc.), concentrated colloidal 



gel, taken from Refs. [3l| and [30], respectively. For 



all systems, data are taken under single scattering con- 
ditions. The experiments on the "artificial skin" will be 
discussed somehow more in depth, because i) its behav- 
ior is representative of the main features observed for the 
other systems; ii) unlike most colloidal systems where 
the control parameters for the slow dynamics are volume 
fraction and interparticle potential, in "artificial skin" 
the slow dynamics is controlled by the relative humidity 
(RH). In this respect, "artificial skin" is representative 
of a class of materials of biological relevance whose slow 
dynamics are still poorly characterized and whose visco- 
elastic properties change dramatically with changes in 
hydration, as in stratum corneum, the uppermost part of 
the epidermis of mammalian skin |4ll443j . This behavior 
is due to the fact that, for these materials, water acts as 
a plasticizer, lowering their glass transition temperature 
below room temperature at high RH. 

The "artificial skin" is a protein-based thin film (20 — 
25 /im) that mimics the properties of human stratum 
corneum. We mount a piece of the film on a circular 
frame of 10 mm of diameter to keep it flat and place 
it in a sealed custom-made cell that allows the relative 
humidity (RH) to be controlled during the PCI measure- 
ments. The desired value of the RH is imposed by plac- 
ing a saturated solution of a suitable salt in a reservoir 
contained in the cell. We present here results for two sat- 
urated salt solutions (prepared by adding 20% more in 
weight fraction than the solubility limit at 20°C), yield- 
ing a humid (RH — 62%, using KI) and a dry (RH = 
12%, using LiCl) atmosphere, respectively. The mea- 
surements are performed in the backscattering geome- 
try depicted schematically in Fig. lb (where the RH- 
controlling cell is not shown for simplicity), with 6* = 180 
degrees, corresponding to a typical probed length scale 
^ ^ 30 nm. The imaged portion of the sample has a 
size of 4.288 x 4.288 mm^ (3.752 x 3.752 mm^) for the 
dry (humid) sample; in both cases, square ROIs of size 
268 X 268 /im^ were used for the PCI analysis. More de- 
tails on the lags used and the duration of the experiment 
are given in the discussion below. 

The colloidal hard sphere samples are suspensions of 
poly-(methyl methacrylate) spheres of radius ~ 100 [45j 
nm in an index matching mixture of cis/trans-decalin and 
tetralin. Their average dynamics have been reported in 
Refs. [13, El]- For the space- resolved measurements dis- 
cussed here, we use the PCI geometry shown in Fig. [T^, 
where g = 25 /im~^, and where the size of the field of 
view and ROIs are 2.5 x 0.6 mm^ and 55 x 55 /zm^, re- 
spectively. We present data for two volume fractions, 
(fi = 0.5468 (delay time r = 20 ms, relaxation time of the 
dynamic structure factor Tq, = 0.147 s) and ip = 0.5957 
(r = 14 s, Ta = 2350 s). All data are taken in a regime 
were the dynamics are stationary. 

The onion gels are concentrated surfactant solutions 



forming a dense packing of polydisperse, deformable 
spheres of average size of the order of a few /im j46j . 
Each sphere is constituted by a stacking of surfactant 
bilayers that roll up to form a multilamellar vesicle, or 
onion. The average dynamics and the rheologic al p rop- 
erties of the onions have been described in Refs. [47l. l48| : 
an optical microscopy investigation [4^ [s^l suggests that 
their dynamics is correlated at least up to length scales 
comparable to the field of view accessible in those exper- 
iments, about 1 mm. In the PCI measurements reported 
here, we use the 90 degree scattering angle geometry of 
Fig. [T^, with a field of view of 1.9 x 0.39 mm^ and ROIs 
of size 130 X 130 fj.m^. The delay time is r = 1000 s, 
about 20 times smaller than the relaxation time of 52 — 1- 
The dynamics slow down with time: we analyze data for 
20000 < t < 60000 s {t = being the time at which 
the surfactant solution is quenched in the gel phase by a 
temperature jump 47]). 

Laponite RD (Rockwood, US) is a synthetic clay con- 
sisting of discoid charged particles. By dispersing 3.5%wt 
of Laponite powder into pure water, we obtain a col- 
loidal glassy suspension that keeps aging for several 
days [5l|, |52|. To properly weight the Laponite content, 
we dry the powder for 12 hours in an oven at 130 C°. In 
order to prevent chemical corrosion, the Laponite powder 
is then dispersed in a pH 10 solution constituted of pure 
Millipore water and a proper amount of sodium hydrox- 
ide. The white dispersion is held in a magnetic stirrer for 
25 minutes until it becomes almost transparent, and it is 
then injected into the scattering cell through a 1 /im filter 
in order to eliminate particle clusters. The sample is pre- 
pared in an inert argon atmosphere in order to prevent 
contact with CO2. The scattering cell is then sealed in or- 
der to preserve this pure condition. The final step before 
measurements consists of a 3 minute centrifugation at 
3000 rpm, to eliminate small gas bubbles that may have 
formed during the process. This sets the age zero of the 
sample (t — 0) with an incertitude of a few minutes. The 
data shown here refer to 2.0 x 10^ sec <t <2.5x 10^ sec; 
they were taken using a special cell as described at the 
end of the next section. 



III. RESULTS 

As a representative example of dynamical heterogene- 
ity in jammed systems, we first describe in some detail 
the behavior of "artificial skin" films. We show in Fig. [2] 
both the average dynamics and the dynamical suscepti- 
bility for a film in humid (RH = 62%, a) and dry (RH = 
12%, b) atmosphere. The data shown here are obtained 
by processing the full images and by averaging over sev- 
eral tens of thousands of sec (36000 sec for the sample at 
RH = 12% and 22500 sec for RH = 62%). The temporal 
intensity correlation functions (black squares, left axis) 
exhibit a decay on a very long time scale, of the order 
of 10"* sec. For the dry sample, §2 — ^ appears to be 
more noisy than for the humid one, with an additional 
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FIG. 2: Intensity time autocorrelation function, (?2 — 1 (left 
axis and solid black squares), and dynamical susceptibility, x 
(right axis and open circles), for the "artificial skin" sample, 
as measured in the backscattering configuration schematized 
in Fig. [TJj. a: relative humidity = 62%; b: relative humidity 
= 12%. Note the enhanced dynamical heterogeneity in dry 
atmosphere, as shown by the higher values of x- The three 
lags indicated by open stars are those for which the spatial 
correlation of the dynamics is shown in FigU) 



decay at earlier times, r « 400 sec. This "noisiness" is 
due to enhanced dynamical heterogeneity, as compared 
to the humid sample: as discussed below, the dynamics 
of the dry sample are intermittent in time and corre- 
lated in space over very large distances: as a result, the 
temporal and spatial averages performed in calculating 
(72 — 1 are poorer than for the humid sample. In both 
panels, we show also the dynamical susceptibility %(''") 
(open circles and right axis). For the humid sample, x 
depends weakly on r, barely exhibiting a peak on the 
time scale of the decay of 32 — 1 • A peak i n y is a distinc- 
tive feature of dynamical heterogeneity [ij, [13, IH, H^] , 
while the noise contribution to the dynamical suscepti- 
bility is essentially flat for time delays smaller than the 
decay time of the intensity correlation function [s^, HI] . 
Thus, the data in Fig. [5^ suggest that the dynamics for 
the humid sample are mildly heterogeneous and that the 
noise contribution dominates the measured Xj with the 
exception of the small peak around r = lO'' sec. By 
contrast, the dynamical susceptibility of the dry sample 
presents a well-developed peak, about a factor of three 
higher than the short-r value of the dynamical suscep- 
tibility and more than a decade higher than x at large 
delay times. This suggests that the dynamics of the dry 
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FIG. 3: For a representative ROI, degree of correlation ci 
at a fixed delay r as a function of time and its PDF for the 
"artificial skin" in humid ( RH = 62%, r = 150 sec, panels 
a and b) and dry (RH = 12%, r = 90 sec, panels c and d) 
atmosphere. Note that the dry sample exhibits large temporal 
fiuctuations. 



sample are highly heterogeneous. 

In order to gain further insight on the nature of the 
dynamical heterogeneity of "artificial skin" , we plot in 
Fig. [3] the temporal evolution of the degree of correla- 
tion, c/, calculated for a time delay much smaller than 
the relaxation time of 52 {t — 150 and 90 sec for the 
humid and dry sample, respectively). We focus on small 
delays because this allows us to capture better the im- 
pact of individual rearrangement events, since at larger 
r several events may occur in between two images. The 
signal shown here have been calculated for an individ- 
ual ROI, but are representative of the general behavior 
of all ROIs. For the humid sample, c/ exhibits modest 
fiuctuations symmetrically distributed around its mean 
value (panel a). Moreover, the probability distribution 
function (PDF) of c/, shown in panel b, is close to Gaus- 
sian. We have shown in previous work [2^ [s^, HI] that 
this is typical of a signal dominated by the measurement 
noise, which is particularly high when the degree of cor- 
relation is calculated over the limited number of pixels 
of a small portion of the image, as required by a space- 
resolved analysis. Panels c and d show c/ and its PDF, 
respectively, for the dry sample. In contrast to the hu- 
mid film, sudden drops of the degree of correlation can be 
easily observed, corresponding to sudden rearrangement 
events in the sample. As a result the PDF of the degree 
of correlation is skewed, with a tail associated with the 
intermittent drops of c/. In other systems, the shape of 
the PDF has been shown to be well fitted by a Gum- 
bel distribution [s^, often observed in the statistics of 
rare events in systems with extended spatial and tempo- 
ral correlations 53|. Here, the dynamics were too slow 
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to allow us to accumulate enough statistics to make any 
quantitative statement on the shape of the PDF. We note 
however that the analysis of the temporal evolution of c/ 
and its PDF confirms the differences between the two 
samples discussed in relation to Fig. [5] the dynamics of 
the dry sample are much more heterogeneous than those 
of the humid one. 

Since the dynamics measured on the full field of view 
(Fig. [2]) exhibit significant fiuctuations, we expect them 
to be spatially correlated over macroscopic distances. 
The spatial correlation of the dynamics calculated ac- 
cording to Eqs. ([5]) and ([6]) is shown in Fig. SI for the 
two samples, and for the three time lags indicated by 
open stars in Fig. [2] Quite generally, G4 is significantly 
larger than zero up to distances of the order of 1 mm 
(for the humid sample) or even several mm (for the dry 
sample). This is remarkable, since all relevant structural 
length scales are much smaller. Furthermore, the range 
of the spatial correlations of the dynamics is much longer 
for the dry sample than for the humid one; together with 
the enhanced intermittency under dry conditions seen 
in Fig. |3l this explains why the dynamical susceptibility 
of the dry sample is higher than that of the wet one, as 
shown in Fig. [2j Interestingly, the time lag dependence of 
G4 is different for the two values of RH. At high RH, the 
range of G4 is maximum at intermediate time delays and 
decreases at small and large r. This is similar to what re- 
ported in numerical work on supercooled molecular glass 
formers (40| and experiments on grains suggesting 
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that dynamical heterogeneity progressively build up with 
time, starting from individual events that are relatively 
limited in size. Eventually, at very large time lags, one 
expects G4 to decay on short length scales, since many 
uncorrelated events will have occurred everywhere in the 
sample. The shorter range of GasXt — 9000 sec observed 
in Fig. 2^ probably corresponds to the onset of this large 
T regime. For the dry sample in panel b, Gi is highly cor- 
related over a very long range for essentially all probed 
lags, indicating that individual events extend over several 
millimeters. We don't have currently an explanation for 
the somehow faster decay of G4 at intermediate lags; the 
long range of 04 at the largest probed lag (t = 9000 sec) 
indicates that on this time scale the regime where many 
uncorrelated events have occurred throughout the sam- 
ple is not yet attained. Longer lags could not be reliably 
analyzed due to the limited duration of the experiment. 

Summarizing, the dynamics of the "artificial skin" 
are markedly heterogeneous under dry conditions, asso- 
ciated with intermittent rearrangement events. These 
events imply the motion of the scatterers over just a few 
nanometers (we recall that the full decay of 52 — 1 corre- 
sponds to displacements of the order of 30 nm); however, 
whenever they occur, they impact essentially the whole 
sample, since the spatial correlation of the dynamics ex- 
tends over several mm. For the humid sample, dynamical 
heterogeneities are less pronounced; the range of spatial 
correlations, although reduced in comparison to that of 
the dry sample, is still on the order of 1-2 mm, much 
larger than any structural length scale in the sample. The 
change in dynamical behavior under different RH condi- 
tions observed here for the "artificial skin" is related to 
the loss of rigidity upon hydration and is reminiscent of 
the sensitivity to humidity of biological tissues. Indeed, 
humidity plays a key role in controlling biomechanical 
properties of stratum corneum [54], whose elastic mod- 
ulus, e.g., increases by several orders of magnitude with 
decreasing RH [H, . 

One may wonder whether similar dynamical features 
are observed also in other glassy or jammed soft matter 
systems that, similarly to the "artificial skin" , are charac- 
terized by ultraslow dynamics and a viscoelastic behavior 
where the solid-like response dominates. To address this 
question, we plot in Fig. [S] the spatial correlation func- 
tion of the dynamics of a series of systems: concentrated 
hard spheres at the onset of the supercooled regime and 
in the deep supercooled regime, soft spheres below and 
above random close packing (data taken form Ref. |31|). 
an onion gel, a colloidal gel (data taken from Ref. |30]), 
and the dry "artificial skin" film already shown in Fig. [21 
With the exception of the soft spheres and the colloidal 
gel, for which data were taken in a low-angle configura- 
tion {q ^ 1 /J,m~^), all data are collected either at 6* = 90 
or 180 degrees. For all samples, we calculate the spatial 
correlation of the dynamics at a time lag r much shorter 
than the relaxation time of the corresponding g2 — l^ in 
order to capture as much as possible the characteristics 
of single events. Quite generally, very long range correla- 
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FIG. 5: Spatial correlation of the dynamics for all systems 
that we have studied except for Laponite, as indicated by the 
label (On: onion gel; CG: colloidal gel, from Ref. [s^]; AS: 
"artificial skin", r = 300 sec; SoS: soft spheres, from [sil ]: 
HS: hard spheres). Data for Laponite are shown in Fig. [6] 
below. All spatial correlation functions have been normalized 
according to Eq. except for the hard spheres, for which 
G4 is shown (see Eq. ((5])), rather than G4. 



tions of the dynamics are observed for all samples, up to 
several mm, once again much larger than any structural 
length scale and comparable to the system size. In partic- 
ular, data for the onions, albeit somehow noisy due to the 
relatively small number of available pixels and the lim- 
ited duration of the experiment (about twice the relax- 
ation time of 52 — 1), confirm and extend the long-range 
correlation of the dynamics observed, for a smaller field 
of view, by microscopy [i^, [s^l . Even for the "artificial 
skin" under humid conditions and the less compressed 
soft sphere system (T = 28°C), where Ga decreases more 
rapidly than for the other systems, spatial correlations 
still extend over macroscopic length scales. The only ex- 
ceptions are the hard spheres samples. For these systems, 
Gi drops to zero as soon as Ar > 0; although we show 
here data at r much smaller than the relaxation time of 
(72 — 1, we point out that a similar behavior is found at all 
T. Since the normalization method described in Sec. Ill Al 
can not be applied, we represent G4 rather than G^ as for 
the other samples. For the less concentrated hard sphere 
suspension, such a sharp drop of G4 is not surprising, 
since no spatial correlation of the dynamics are to be ex- 
pected at the onset of the supercooled regime. Indeed, 
at = 0.5468 the dynamics is slower than in the ip ^ 
limit by a relatively modest factor of 100, and the shape 
of the dynamic structure factor differs only marginally 
from that in the diluted case [2Q|. For the most con- 
centrated hard sphere sample, by contrast, the system 
relaxation is already slowed down by a factor of about 
10^ compared to the (p 7^ Q limit and a fully developed 
plateau is observed in the dynamic structure factor, re- 
vealing caging and glassy dynamics [lO, El] . It has been 
proposed that such a dramatic slow down of the dynam- 
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FIG. 6: Bottom left: schematic side view of the cell used 
for the experiment on Laponite. LB: laser beam, GS: glass 
slide. Top: typical CCD image of the scattering volume. The 
dark region corresponds to the thickness of the glass slide, 
view from the side. The three ROIs for which the degree of 
correlation is shown in the bottom right plot are highlighted. 
The size of the imaged region is 2.62 x 0.52mm^. Bottom 
right: time dependence of c/, for a delay time r = 200 sec, 
for the three ROIs shown above. For the sake of clarity, the 
curves of ROIs 1 and 2 have been offset vertically by 0.2 and 
0.1, respectively. Note that the signals measured on the same 
side of the glass slide (ROIs 1 and 2) are correlated, while 
signals from opposite sides are uncorrelated. 



ics is associated with the growth of spatial correlations 
of the dynamics. Previous experiments, both by confo- 
cal microscopy [l^ and dynamic light scattering [^, [lO] , 
indicate that the range, ^, of such correlations does not 
exceed a few particles sizes, i.e. a few fiia in our case. 
Therefore, the absence of any measurable correlation in 
our PCI experiment, where the smallest accessible spatial 
lag is 55 iim, is fully consistent with previous results. On 
the one hand, this negative result illustrates the limita- 
tions intrinsic to PCI, a coarse grained method. On the 
other hand, it highlights the dramatic difference between 
the dynamics of supercooled, stationary hard spheres and 
those of the other jammed, out-of-equilibrium samples, 
for which macroscopic spatial correlations of the dynam- 
ics are observed. 

The results shown here suggest that system-size corre- 
lations of the dynamics are ubiquitous in fully jammed 
systems, where the solid-like behavior of the sample dom- 
inates over its viscous response. Indeed, in a material 
that is essentially solid-like, any deformation due to a 
local rearrangement will propagate very far before being 
substantially damped, leading to extended spatial corre- 
lations of the dynamics. Another (less attractive) expla- 
nation of the dynamical fluctuations reported here could 
be that they are due to some kind of artifact. In fact, 
one should be aware that measuring fluctuations in a re- 



liable fashion is much more difficult than just probing the 
average dynamics, and many spurious effects have to be 
controlled and ruled out (for a discussion of some possi- 
ble sources of artifacts, such as temperature fluctuations 
and laser beam pointing instabilities, see also Ref. ^38,] ) . 
In order to address the role of the elastic propagation of 
a strain fleld and to demonstrate the genuine nature of 
the spatial correlations observed here, we discuss briefly 
an experiment performed on a Laponite sample using a 
specially devised cell. 

The cell is a cylindrical tube of inner diameter « 10 
mm (see Fig. [6|) , with a microscope glass slide glued in 
the interior with epoxy, so as to separate the cell into 
two chambers. Because the glass slide is cut unevenly 
and does not fit perfectly the cell walls, the two cham- 
bers arc in contact through openings of size « 1 mm, 
so that the sample contained in the two chambers is 
fully equilibrated and experiences exactly the same con- 
ditions. The laser beam impinges perpendicularly to the 
glass slide, thereby illuminating both chambers, and the 
scattering volume is imaged at an angle 9 = 90°, as in 
the geometry shown in Fig.[T^. We calculate the local cj 
for larger-than-usual ROIs, in order to reduce the noise 
contribution (see Fig. [6|). From the plot in Fig. \6\ it is 
clear that the time evolution of the degree of correlation 
measured on the same side of the glass slide is highly 
correlated (ROIs 1 and 2, separated by Ar = 0.57 mm), 
since most of the downward spikes corresponding to sud- 
den rearrangement events are co-occurrent. Indeed, G4 
calculated from the signals for ROIs 1 and 2 is as high as 
0.394 (we recall that G4 is affected by the noise contribu- 
tion, so that an even higher value would be observed in 
the absence of noise). By contrast, the events recorded 
in ROI 3, on the opposite side of the glass slide, appear 
to be uncorrelated with respect of those of ROIs 1 and 2, 
as confirmed by low values of the correlation (e.g., = 
-0.015 when correlating the signals from ROIs 2 and 3). 

This experiment demonstrates that the downward 
spikes observed in c/ can not be due to some mechanical, 
temperature or laser beam instability, since they occur 
independently on either side of the glass slide. On the 
other hand, the data of Fig. [S] support the idea that the 
propagation of a strain field in a predominantly elastic 
medium is responsible for long range spatial correlations. 
Indeed, such a strain field can not propagate through the 
rigid glass slide, so that the dynamical activity in the 
two chambers is uncorrelated, while the one in a given 
chamber is correlated. 



IV. DISCUSSION 

To better understand the relationship between vis- 
coelastic properties and spatial correlations of the dy- 
namics, we examine published rheological data for most 
of the systems investigated here and summarize the rele- 
vant rheological parameters in Table [H where we choose, 
somehow arbitrarily, j/ ~ 1 Hz as a reference frequency 
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TABLE I: Rheological parameters of most of the system 
shown in Figs. [S] and [B] The sample names are as in the 
caption of Fig. [S] 





V (Hz) 


G'{v) (Pa) 


G'{^)/G"{v) 


Ref. 


On 


1 


600 


15 


m 


CG 


1 


~ 0.9 X 10~^ 


10 


[59] 


SoS, ^ = 0.57 


1.6 


0.6 


0.3 


m 


SoS, if = 0.69 


1.6 


20 


8 




HS, LP = 0.5468 


1 


40 


1.1 


[60, 61] 


HS, (/p = 0.5957 


1 


> 80 


> 1.4 


[60, 61J 


Laponite 


0.7 


> 300 


20 


[55, 56] 



for comparing different systems. Oscillatory rheology 
experiments for the onion gels [4^ show that both the 
storage, G'^u), and loss, G"{v), moduli are essentially 
frequency-independent in the range 4 x 10"'^ Hz < v < 
10 Hz, and that G' ^ 600 Pa dominates over G" by a fac- 
tor of ^ 15. Similar results are obtained for the Laponite 
suspension and for the colloidal gels [F^. Note that 
the magnitude of the elastic modulus of the colloidal gel 
(G" - 4.4 X 10--* - 1.5 X 10-3 Pa is lower than 

that of the onion gel and of the Laponite suspension by at 
least 5 orders of magnitude, implying that the absolute 
value of G' is not a relevant parameter in determining the 
range of spatial correlations of the dynamics. A system- 
atic investigation of the shear moduli of the soft particles 
is presented in Ref. [sTI . For the most concentrated sam- 
ple reported in Fig. [S] above, one finds that G" and G" 
are nearly frequency independent and that the former is 
larger than the latter by almost a decade. By contrast, 
for the less concentrated sample both moduli are com- 
parable (see panels a) and b) of Fig. 3 of Ref. [3l|). 
revealing a complex viscoelastic behavior where neither 
the solid-like nor the fluid-like character prevail. Since 
spatial correlations of the dynamics are shorter-ranged 
in the less concentrated sample, these observations are 
consistent with the notion that a fully developed elastic 
behavior is a necessary condition for long-range spatial 
correlations of the dynamics. The same trend is likely 
to hold for the "artificial skin". Indeed, we find that 
the range of the spatial correlation of the dynamics in- 
creases as the RH decreases. Although no rheology data 
are available for the "artificial skin", it is worth noting 
that in stratum corneum [4^, |4^ , a material whose prop- 
erties the Vitro-corneum® film is designed to mimic, the 
elastic modulus grows by almost 3 orders of magnitude 
when the RH decreases from 100% to 30%. 

Collectively, these results indicate that systems with 
long-ranged spatial correlations of the dynamics have a 
solid-like behavior. A close inspection of rheology data 
for hard sphere suspensions, for which we recall that ^ 
is limited to a few particle sizes at most, show that the 
reverse is not true. Mason and Weitz have measured 
the shear moduli of supercooled hard spheres (60j . Care 
must be taken in applying their results to our hard sphere 
samples, due to the unavoidable uncertainties in the de- 



termination of the absolute volume fraction [3, . A 
generally accepted approach consists in considering the 
separation, e — [kpc — ip)/ipc, with respect to the location 
Lfc of the (apparent) divergence of the relaxation time as 
obtained from a mode coupling theory fit of (i p) [ In our 
work, w 0.59 H, while ~ 0.575 in Ref. For 
our most diluted sample, Lp = 0.5468 and e « 0.07, cor- 
responding roughly to the data for = 0.53 in Ref. [60t . 
At that volume fraction, G' and G" have approximately 
equal magnitude in all the accessible frequency range, 
similarly to the case of our soft spheres dX — 0.57. How- 
ever, no spatial correlation of the dynamics is observed 
in the hard spheres aX ip = 0.5468, while ^ 2 mm for 
the soft spheres a,t (p = 0.57. This difference is even more 
striking if we consider our most concentrated hard sphere 
sample, for which ^ is limited to a few particle sizes poj . 
although the sample is solid-like (G' dominates over G") 
and the magnitude of its elastic modulus [6l| is compa- 
rable to that of the Laponite suspension, which exhibits 
systems-size dynamical correlations. 

Our measurements thus indicate that, while elasticity 
is a necessary condition for observing extended spatial 
correlations of the dynamics, it is not a sufficient one. 
The microscopic origin of the elasticity must also play a 
crucial role: the differences discussed above suggest that 
^ is modest in systems where the elasticity has an en- 
tropic origin, as for the hard spheres, and no long-lived 
network of contacts exists. By contrast, ^ is very large 
in systems where elasticity has an enthalpic origin, be it 
due to the bending rigidity of the protein filament net- 
work in "artificial skin" and the backbone in diluted gels 
made of strongly attractive particles, or to the bulk elas- 
ticity of squeezed particles, as in the soft spheres and the 
onions. Note that the Laponite suspension falls also in 
the latter category: while the platelets are not in direct 
contact, they interact via a screened Coulombic repulsive 
potential and can thus be effectively regarded as squeezed 
soft particles. Interestingly, recent work [6^ has identi- 
fied a structural length scale that controls the mechani- 
cal properties of attractive colloidal systems: this length 
scale is of the order of the particle size for hard spheres, 
while it grows as short-ranged attractive interactions be- 
come increasingly important, as in colloidal gels. Our 
data suggest that this behavior may be mirrored by a 
similar growth of ^. In repulsive, squeezed systems, the 
same role is presumably played by a persistent network 
of interparticle forces, such as that visualized by confocal 
microscopy in compressed emulsions f63 |. 

In view of the above discussion, it is natural to com- 
pare our results to numerical and experimental work on 
driven jammed systems. In these works, the sample is 
sheared by imposing either a continuous deformation or 
a sinusoidal one. The dynamics are quantified by the par- 
ticle displacements perpendicular to the shear direction 
or after subtracting the affine component (for a continu- 
ous shear), or by comparing successive configurations at 
zero deformation (for oscillatory shear). For a 2D gran- 
ular medium sheared at a finite shear rate, Lechenault 
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et al. [34I find that a spatial correlation function analo- 
gous to our 6*4 decays on a length scale ^4 on the order 
of ten particle sizes at the jamming transition, and that 
this length decreases upon further compression, although 
only a restrained range of densities above jamming could 
be probed, due to the particle stiffness. Simulations of 
sheared soft disks 333 indicate that above the jamming 
transition and at finite shear stress, a, ^ remains moder- 
ate, but a scaling analysis suggests that ^ diverges in the 
limit (7 — > for all densities above jamming. This is con- 
firmed by quasi-static shear simulations of jammed soft 
particles [sj, , where it is shown that spatial correla- 
tions of the dynamics above jamming are limited only by 
the system size. 

Therefore, the dynamical behavior reported here ap- 
pears to be due to the connected nature of the materials 
investigated, where a strain field can propagate over very 
large distances, as in fully jammed systems sheared at a 
vanishingly small rate. A related issue, still open, con- 
cerns the microscopic ori gin of the dynamics. The sys- 
tems studied in Refs. |33 - l35| are athermal, with no dy- 
namics in the absence of an external driving. Although 
the systems presented here are thermal and no external 
drive is applied to them, it is unlikely that thermal mo- 
tion alone is responsible for their dynamics in the jammed 
state. This is illustrated, e.g., by the volume fraction de- 
pendence of the dynamics reported for the soft spheres 
of Ref. [3l|. In the supercooled regime, the relaxation 
time, Tc, grows sharply with ip, as in hard sphere systems. 
Above random close packing, a different regime sets in, 
where Tc grows very slowly with if. Thus, above jamming, 
the relaxation time is orders of magnitude smaller than 
what expected by extrapolating the behavior in the su- 
percooled regime, strongly suggesting that an additional 
relaxation mechanism has set in. We propose that, quite 
generally, the relaxation of internal stress may be such 
a mechanism in jammed soft matter. This is an appeal- 
ing explanation, since it would be consistent with the 
observed ultra-long range correlations of the dynamics. 
However, this conjecture still awaits for a direct experi- 
mental proof. 

V. CONCLUSIONS 

We have presented direct measurements of the corre- 
lation length of the slow dynamics of a variety of glassy 



and jammed soft systems, obtained using the recently 
introduced PCI method. This technique allows one to 
measure coarse grained maps of the dynamical activity 
of a sample. Its main advantage consists in the possibil- 
ity of probing motion on a very small length scale, yet 
for a very large field of view. This is an attractive feature 
for jammed systems, where motion is very restrained but 
highly correlated spatially. Additionally, the technique 
does not require individual particles to be imaged; ac- 
cordingly, it can be applied to systems for which direct 
visualization by optical or confocal microscopy is not pos- 
sible, including turbid samples as, e.g., foams [65| . 

We find that in deeply jammed systems ^ is quite gen- 
erally very large, typically on the order of the system 
size. This is in striking contrast with the modest corre- 
lation lengths measured in glass formers, including col- 
loidal hard spheres. An analysis of the viscoelastic prop- 
erties of the various systems investigated here shows that 
a well developed elasticity (C > G") always accompa- 
nies the presence of long-ranged spatial correlations of 
the dynamics, regardless of the absolute magnitude of 
the elastic modulus. The reverse, however, is not true, 
as exemplified by concentrated hard sphere suspensions, 
whose macroscopic rheological response is predominantly 
solid-like, whereas the correlation length of the dynamics 
is limited to a few particle sizes. These results highlight 
the crucial role of the microscopic origin of elasticity (en- 
tropic vs. enthalpic) in determining the range of spa- 
tial correlations of the dynamics. Further work will be 
needed to test thoroughly these ideas, for example by 
varying continuously the strength of attractions in dense 
colloidal suspensions, so as to change progressively the 
nature of the elasticity [6^ and, presumably, the corre- 
lation length ^. 
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